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Abstract 

We study theoretically the space-time evolution of the thermal and electromagnetic perturbation in a 
superconductor with a nonlinear current-voltage characteristics in the flux creep regime. On the basis 
of a linear analysis of a set of differential equations describing small perturbations of temperature 
and electromagnetic field, it is found that under some conditions a branching instability may occur 
in a superconductor sample. 
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The dynamics of thermomagnetic instabilities of the critical 
state and flux jumps in hard superconductors with high values of 
the critical current density and the critical magnetic fields have 
been earlier investigated by many authors [1-5]. The general con- 
cept of thermomagnetic instabilities in type-II superconductors 
was developed in literature [5]. The dynamics of small thermal 
and electromagnetic perturbations, whose development leads to 
the flux jump, have been investigated theoretically by Mints and 
Rakhmanov [4, 5]. The authors have found the stability criterion 
for the flux jumps in the framework of adiabatic and dynamic 
approximations in the viscous flux flow regime of type-II super- 
conductors. The detailed theoretical analyze of the flux jumping 
in the flux creep regime, where the current-voltage characteris- 
tics of a sample is a nonlinear have been carried out recently by 
Mints [6] and by Mints and Brandt [7]. Nonuniform magnetic 
flux penetration in type II superconductors, creating finger and 
dendritic patterns, has recently attracted considerable interest. 
Such patterns have been directly observed in a large number of 
superconducting films employing magneto-optical imaging tech- 
niques [8-15]. The existing experimental data [8-15], and the 
recently developed theoretical models [16, 17], suggest that the 
origin of these patterns is thermomagnetic instability of the vor- 
tex matter in the superconducting films. The instability arises 
from local temperature increase due to flux motion, which, in 
turn, decreases flux pinning and hence facilitates further flux mo- 
tion. Linear stability analysis, based on the coupled nonlinear 
Maxwell and thermal diffusion equations, showed that instability 
in the form of narrow fingers perpendicular to the background 
electric field occurs when this field exceeds of its threshold value 
[16, 17]. 

In the present work, we study the spatial and temporal evolu- 
tion of small thermal and electromagnetic perturbation in type-II 
superconductor sample in the flux creep regime with a nonlinear 
voltage-current characteristics, assuming that an applied field 
parallel to the surface of the sample. On the basis of a linear 
analysis of a set of differential equations describing small pertur- 
bations of temperature and electromagnetic field we found that 
under some conditions a branching instability may occur in the 
sample. 

Mathematical problem of theoretical study the dynamics of 
thermal and electromagnetic perturbations in a superconductor 
sample in the flux creep regime can be formulated on the basis of 
a system nonlinear diffusion-like equations for the thermal and 
electromagnetic field perturbations with account nonlinear rela- 
tionship between the field and current in superconductor sample. 
The distribution of the magnetic flux density B and the trans- 



port current density j inside a superconductor is described by 
the equation 

rotB = no j. (1) 

When the penetrated magnetic flux changes with time, an elec- 
tric field E is generated inside the sample according to Faraday's 
law 

rot E = ^L. (2) 
dt y ' 

The temperature distribution in superconductor is governed by 

the heat conduction diffusion equation 

v(T) d ^=V[K(T)VT]+jE, (3) 

Here v — v(T) and k = k(T) are the specific heat and ther- 
mal conductivity, respectively. The above equations should be 
supplemented by a current-voltage characteristics of supercon- 
ductors, which has the form 

j=j c (T,B)+j(E). 

In order to obtain analytical results of a set Eqs. (l)-(3), we sug- 
gest that j c is independent on magnetic field induction B and use 
the Bean critical state model j c — j c (B e , T) = jo — a(T — To) [1, 
18], where B e is the external applied magnetic field induction; 
a = Jo/(T c — To); jo is the equilibrium current density, To and T c 
are the equilibrium and critical temperatures of the sample, re- 
spectively, [5]. For the sake of simplifying of the calculations, we 
perform our calculations on the assumption of negligibly small 
heating (T — To <§C T c — To) and assume that the temperature 
profile is a constant within the across sample and thermal con- 
ductivity k and heat capacity v are independent on the temper- 
ature profile. We shall study the problem in the framework of a 
macroscopic approach, in which all lengths scales are larger than 
the flux-line spacing; thus, the superconductor is considered as 
an uniform medium. The system of differential equations (l)-(3) 
should be supplemented by a current- voltage curve j = j(E). In 
the flux creep regime the current- voltage characteristics of type 
II conventional superconductors is highly nonlinear due to ther- 
mally activated dissipative flux motion [19, 20]. For the logarith- 
mic current dependence of the potential barrier U(j), proposed 
by [21] the dependence j(E) has the form 



where Eo is the voltage criterion at which the critical current 
density j c is determined [5]; a constant parameter n depends 



1 



on the pinning regimes and can vary widely for various types of 
superconductors. In the case n = 1 the power-law relation (4) 
reduces to Ohm's law, describing the normal or flux-flow regime 
[18]. For infinitely large n, the equation describes the Bean crit- 
ical state model j = j c [1]. When 1 < n < oo, the equation (4) 
describes nonlinear flux creep [22]. In this case the differential 
conductivity a is determined by the following expression 



(5) 



a = ^- = 

dE nE b 

According to relation (5) the differential conductivity decreases 
with the increasing of the background electric field Eb, and 
strongly depends on the external magnetic field sweep rate 
Eb ~ B e x. Therefore the stability criterion also strongly de- 
pends on the differential conductivity a. For the typical values 
of ji = 10 3 A/cm 2 , E b = lQ- 7 V/cm we obtain a = 10 u) l/Q.cm. 
It follows from this estimation [6, 7] that the differential con- 
ductivity a which determines the dynamics of thermomagnetic 
instability for is high enough. We assume, for simplicity, that 
the value of n temperature and magnetic-field independent. 

Let us formulate a differential equations governing the dy- 
namics of small temperature and electromagnetic field pertur- 
bation in a superconductor sample. We study the evolution of 
the thermal and electromagnetic penetration process in a simple 
geometry - superconducting semi-infinitive sample x > 0. We 
assume that the external magnetic field induction B e is parallel 
to the z-axis and the magnetic field sweep rate B e is constant. 
When the magnetic field with the flux density B e is applied in the 
direction of the z-axis, the transport current j(x, t) and the elec- 
tric field E(x, t) are induced inside the slab along the y-axis. For 
this geometry the spatial and temporal evolution of small ther- 
mal T(x, t) and electromagnetic field E(x, t) perturbations are 
described by the thermal diffusion equation coupled to Maxwell's 
equations 



' dt 



d 2 E 
dx 2 



jc dE 
nEb dt 



dj c dT 
dT~dt 



(6) 



(7) 



It should be noted that the nonlinear diffusion-type equations 
(6) and (7), totally determine the problem of the space-time dis- 
tribution of the temperature and electromagnetic field profiles in 
the flux creep regime with a nonlinear current-voltage character- 
istics in the semi-infinite sample. 

Let us specify the thermal and electrodynamic boundary and 
initial conditions to the last system of equations. The thermal 
boundary conditions are 



dT{Q,i) 
dx 



0, T(L,t) = T 



We assume that the magnetic field perturbation is equal to zero 
at the sample surface and according to relation (2), we obtain 
the first electrodynamic boundary condition 



dE(0,t) 
dx 



0. 



The second boundary condition for the electric field E(x, t) at 
the flux front x = L can be presented as 

E(L,t)=0, 



The boundary conditions for the magnetic induction are 
dB(0,t) = B e , B(L,t) = 0, 

cB e 



where L 



is the London penetration depth. For initial 



conditions we assume that the electric field is uniform within the 
cross-section of the sample E = Eo at t=0. 

Let us derive, for this geometry a differential equations, de- 
scribing the spatial and temporal evolution of thermal T(x,t) 
and electromagnetic field E(x, t) perturbations. We present the 
small thermal and electromagnetic perturbations in the form 



Q(x, t) = T (x) + (T c - To)6 exp [yt/to + iqz] , 
e(x,t) = E b (x) + E b eexp [yt/to + iqz] . 



(8) 



where To(x) and Eb(x) are solutions to the unperturbed equa- 
tions obtained in the quasi-stationary approximation describing 
the background distributions of temperature and electric field in 
the sample. Here 7 is the eigenvalue of the problem to be deter- 
mined and q — 2nd/ L is the wave-number of the perturbation. 
From solutions (8), one can see that the characteristic time of 
thermal and electromagnetic perturbations t is of the order of 
to/7. Where, we have introduced the following dimensionless 
parameters and variables 



to 



ava x , va it d 

— , z =- v d = -, 1=oY- 

Jc d nojc 2 L 

As we mentioned above, the background temperature To(x) is 
practically uniform over the cross-section of the sample and un- 
der this approximation we ignore its coordinate dependence. We 
note that the background electric field may be created by ramp- 
ing the external magnetic field, and for simplicity we assume it 
to be coordinate independent. It turns out that these simplifi- 
cations have no qualitative influence on the results but make it 
possible to perform analytical calculations completely. 

Substituting the expression (8) into the system equations (6), 
(7) one can get the following linearized system equations for O 
and e 



T<2 2 e + 7 e+ -@-2 ( - } e = 0, 
n \ n , 



q 2 e + 7 [e - n6] = 0. 



(9) 



Solving the above system equations (9) we obtain the following 
dispersion relation to determine an eigenvalues of the problem 



7 + 



(r + l)q 2 



7 + 



2 , 
rq + 



= 



(10) 



where r is the ratio between the characteristic time of magnetic 
flux diffusion and the characteristic time of heat flux diffusion 
[4]. The instability of the flux front is defined by the positive 
values of the growth rate Re 7>0. It can be seen that there is a 
critical wave number, 



9c 



(11) 



below which the system is always unstable at n=l. This insta- 
bility appears first at q=0. In this case the small perturbations 
grow with the maximal possible rate Re 7 = 1. The growth rate 
dependence on the wave number for different values of r is il- 
lustrated in Figs. (1-3) at n = 1. For high enough values of 



2 



t the system is stable. As the r decreases, the growth rate 7 
increases. The branching instability will gradually appears for 
relatively small values of t=0.05. 

Thus, according to (11), the branching instability occurs at 
the threshold electric field E = E c 



E = E C = 



ty 2 k(T c - To) 



4 j c L 2 • 

Taking into account an expression for penetration depth L, the 
threshold field can be written at B — B t h as 



(12) 



7T K(T C -T )j c 

Bth= 2^ E b • 

The threshold field for branching instability, as can be seen from 
the last expression is highly sensitive to the critical current den- 
sity and the shape of the background electric field Eb generated 
by the varying magnetic field. The threshold field B t h decreases 
monotonously with increasing the background electric field Eb. 
If we assume that the background electric field Eb generated by 
a varying magnetic field as E b ~ B e , for the considered simple 
geometry, then we can easily obtain the expression for the sweep 
rate dependence of the threshold field of branching instability. 

Let us assume that the thermal diffusion is slower than the 
magnetic diffusion r <C 1. In this adiabatic limiting case, the 
instability criterion is determined as q = q c — 1 for n=l, so the 
threshold field can be presented as 



(Tc 



This is a well known adiabatic criteria [2], which assumes that 
the heat transport from the sample surface to the environment 
can be neglected. 

Conclusion 

In conclusion, on the basis of a linear analysis of a set of differ- 
ential equations describing small perturbations of temperature 
and electromagnetic field we found that under some conditions 
a branching instability may occur in the sample. 
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Fig. 1-3. The dependence of the growth rate on the wave number for r = 0.5, 0.05, 0.08 and n=l. 
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